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Analytical treatment of skyrmions given by rational map (RM) ansaetze proposed recently for the 
Skyrme model is extended for the model including the 6-th order term in chiral field derivatives 
in the lagrangian and used for the calculations of different properties of multiskyrmions. At large 
baryon numbers the approximate solutions obtained are similar to the domain wall, or to spherical 
bubbles with energy and baryon number density concentrated at their boundary. Rigorous upper 
. bound is obtained for the masses of RM multiskyrmions which is close to the known masses, 

' especially at large B. For the 6-th order variant the lower bound for masses of RM skyrmions is 

^Sj I obtained as well. The main properties of the bubbles of matter are obtained for arbitrary number 

of flavours. They are qualitatively the same for the 4-th and 6-th order terms present in the 
O ' lagrangian, although differ in some details. 

o 

'. 1 Introduction 

(N 
> 



Soliton models of different kinds are used in various fields of physics. In elementary parti- 
es) ■ cle physics the soliton models provide a concept of baryons as extended in space objects, 



O 



opposite to the concept of the point-like objects, usual for the quantum field theory. The 



O '• chiral soliton approach, starting with several basic principles and ingredients incorporated 



in the model lagrangian ^ provides realistic and even satisfactory description of baryons 
and baryonic systems. The latter are obtained within this approach as quantized solitonic 
p i' solutions of equations of motion, characterized by the so called winding number or topolog- 
p • ical charge which is identified with the baryon number B. Numerical studies have shown 
^ ■ that the chiral field configurations of lowest energy possess different topological properties 
- the shape of the mass and B-number distribution - for different values of B. It is a sphere 
for B = I hedgehog a torus for 5 = 2, tetrahedron for 5 = 3, cube for 5 = 4, and higher 
^ i polyhedrons for greater baryon numbers. The symmetries of various configurations for B up 
I to 22 and their masses have been determined in (the references to earlier original papers 
where the symmetries of configurations with smaller baryon numbers have been determined 
can be found in 0, These configurations have one-shell structure and for 5 > 6 all of 
them, except two formed from 12 pentagons and 25-14 hexagons; in carbon chem- 

istry similar structures are known as fuUerenes [^. The mass and baryon number densities 
for these configurations are concentrated along the edges of polyhedrons. All these config- 
urations can be made of 25 - 2 slightly deformed torus-like configurations glued together, 
which can be considered by this reason as elementary building blocks for multiskyrmions. 
As will be shown here, the dimensions of these elementary cells do not depend on 5 when 
5 is large enough. 

The so called rational map (RM) ansatz, proposed for the SU{2) skyrmions in 
and widely used now, in present paper as well, allows to simplify the problem of finding 
the configurations of lowest energy. For the RM ansatz the minimization of the skyrmions 
energy functional proceeds in two steps: at first step the map from is minimized 

for the SU{2) model (for the SU{N) model it is a map from S*^ CP^^^, 0), and, second, 
the energy functional depending on skyrmion profile as a function of distance from center of 
skyrmion is minimized. As will be shown here, just the second step can be done analytically 
with quite good accuracy. Many important properties of the RM multiskyrmions can be 



studied in this way, and some of them do not depend on result of the first step. This aUows 
to make certain conclusions for the arbitrary large B and for any number of flavours Np ^ N 
independently of presence of numerical calculations. Without difficulties the consideration 
has been extended to the variant of the model with the 6-th order terms in chiral derivatives 
included into lagrangian (the SK6 variant of the Skyrme model). Remarkably, that for the 
SK6 variant of the model the dependence of the results on the first step of calculation is 
even weaker than for the SK4 variant. 

Beginning with [|l], |^ , the chiral soliton models have been considered as a special class 
of models for baryons. Their connection with other models could be instructive and useful, 
and this is also an issue of present paper. In particular, it is shown that Skyrme-type mod- 
els provide field theoretical realization of the bag model of special kind for baryonic systems. 



2 Large B multiskyrmions as spherical bubbles or domain walls 

Here we consider the multiskyrmions in the general SU{N) case; detailed comparison of 
analytical results with numerical calculations is made in the SU{2) model and also in the 
SU{?,) variant using the projector ansatz 0. In the SU{2) model the chiral fields are functions 
of the profile / and the unit vector n, according to definition of the unitary matrix U e SU{2) 
U = Cf + isfUT. For the ansatz based on the rational maps the profile / depends only on 
the variable r, and the components of vector n - on angular variables 6, 0. ni = {2ReR)/{l + 
I^P), "-2 = (2/TOi?)/(l + na = (1 - |i?P)/(l + where i? is a rational function of variable 
z = tg{9/2)exp{i4>) defining the map of degree TV from 5^. 
The notations are used IIHI 



2idRdR 



4 



J_ f 4[5man2p J_ f (l + lzp) \dR\ y 2idzdz 

where is a spherical angle. For 5 = 1 hedgehog 7V = X= l. N = B for configurations of 
lowest energy. 

For more general SU{N) case and using projector ansatz one obtains 

I dzdzTr[d,Pd-,P], 
2tt J 

J = ^ y dzdz{l + \z\^fTr[d,P, d,P]\ (2) 

P is a projector, iVxA^ hermitian matrix, P = ffyj^f, TrP = 1. For SU{2) case the 2- 
component column / = (i?, 1)^, and (1) can be obtained easily from formulas (2). 

The classical mass of skyrmion for RM ansatz in universal units 'in'^F.j^/e is 0: 

^'^^hj l^^^'-^" + + (1 - A)/'2] + (1 - A)X^ + AX^/'^jdr, (3) 

r measured in units 2/(^Ve), the coefficient Am = 2{N ~ l)/N for symmetry group SU{N) 0, 
and this generalization provides a possibility to consider models with arbitrary number of 



flavours N ^ Np - essentially nonembeddings of SU{2) in SU{N). A defines the weight of 
the 6-th order term. If A = we obtain the original Skyrme model variant (we call it the 
SK4- variant), A = 1 corresponds to the pure SK6- variant. The expression (3) for mixed 
variant was obtained in 0, where the properties of multiskyrmions with B < 5 have been 
investigated in the modification of the Skyrme model with 6-th order term included into 
lagrangian. For 5 = 2 the results obtained in agree with those obtained in which 
correspond to the value of A = 1. Besides axially symmetrical configurations with B < 5 in 
the SK4 variant the torus-like configuration with 5 = 2 was obtained in also in the SK6 
variant. For reasonable values of model parameters the binding energy was about twice 
greater than that for pure Skyrme model: it was found ~ 150 Mew instead of 74 Mew for the 
variant of the SK4 model which provides fit for the masses of nucleon and A-isobar. 

It was important in deriving (3) that the integral over angular variables for the trace 
of the third power of [dzP,dzP] equals to zero 0: 



M = ^ [ dzdz{l + \z\YTr[d,Pd-,Pf = 0, 

OTT J 



therefore, the structure (3) is specific for the rational map ansatz. 

To find the minimal energy configuration at fixed N = B one minimizes J, and then 
finds the profile /(r) by minimizing energy (3). The inequality takes place: I > B^ [H? i|- The 
important consequence of (3) is that the symmetries of multiskyrmions in the SK6 variant 
are the same as for the SK4 variant, because the quantity I is the same for both variants. 

Numerical calculations peformed in 0-0 have shown, and the analytical treatment 
of and here supports, that at large B and, hence, large I the multiskyrmion looks like 
a spherical bubble with profile equal to / = tt inside and / = outside. The energy and 
B-number density of this configuration is concentrated at its shell, similar to the domain 



walls system considered in [|T0[ in connection with cosmological problems. 

The lower bound for the mass of solitons can be obtained from (3) for the SK6 variant, 
similar to the SK4 variant, known previously. Using evident relations A^f^r"^ +Isjf'^/r^ > 
2y/A^f^s), and f'^^/A^ + B > -2f' {B^/A^^'^ (recaU that /' < O) we obtain from (3) for 
arbitrary value of real positive parameter A < 1: 



M{B,X) f AnIY'* ^ 1-A 



B - 3 \ B^ 



1/4 

(4a) 



For A = the known results of ^ are reproduced, for the SK6 variant (A = l) it is new 
bound. When the last term in (3) can be neglected, another lower bound can be obtained, 
see also 0: 



M{B,\) ^ Vl~A 



B - 3 

which does not provide real restriction when A is close to 1. Recall that the bound M/B > 1 
was obtained for arbitrary, not only RM, SU{2) skyrmions first by Skyrme in the SK4 variant 
of the model [|1|. 

Denote = cos/, then the energy (3) can be presented as 

M = ^ II [A^r'^'^ + 2B{1~ + (1 - A) [2B^" + 1{1 <l>r/r'] + 

+XI{l-(j,^)cb'yrAdr, (5) 



with (j) changing from -1 at r = to 1 at r ^ oo. The first part of (5) is the second order term 
contribution into the mass, the second - the Skyrme term contribution, and the last, propor- 
tional to A, - the 6-th order term. At fixed r = ro the 4-th order term is proportional exactly 
to 1-dimensional domain wall energy. It is possible to write the second order contribution 
in (5) in the form: 



Sir J i(l-<^2) 



An 



+ 2r^/2ANB(t)'^dr, 



and similar for the 4-th order Skyrme term. The equality (f)' = ^2B/An{1 - <P^)/r eliminates 
considerable part of integrand for M^'^\ Therefore, it is natural to consider function 4> 
satisfying the following differential equation: 

<^' = A(l-.^2) (6) 

with constant power &, which has solution satisfying boundary conditions (/)(0) = -1 and 
(/>(oo) = 1: 

{r/r^f + 1 

with arbitrary rg - the distance from the origin of the point where and profile / = 7r/2. 
tq can be considered as a radius of multiskyrmion, both b and ro will be defined further by 
means of the mass minimization procedure. The radii of distributions of baryon number 
and mass in the soliton are close to ro, see Section 5 below. 

After substitution of this ansatz one obtains the soliton mass in the form: 



{r.r.M='-;^;^ (7) 



M{B, b) = -^ y"|(Ajv6V4 + 2S)(1 



./>2) + (1 - X){Bb^ /2+I){l - cl>''f/r^+ 



Integrating over dr can be made using known expressions for the Euler-type integrals, e.g. 

io l + {r/rof bsin{n/by ' 

and, more generally 

r WTT7% = ^''"'-'" \ ^ MbV P>o,b>i + c,c>-i. (9) 

Jo 13+ {r/rof 6sin[7r(l + c)/b] 

Differentiation in /? allows to get the integral with any power of 1 -i- {r/rof in denominator. 
What we need now is {(j) is defined in (7)) 



b'^sin{-K/b)^ J 3rob'^sin{'iT/b)' 



(1 - <^2)3 ^ 327r(l-9/6^)(l-9/46^) ^ 32n ^^^^^^ ^^^^ 

brgb"^ sin{'iTT / b) 15rQb'^sin(Tr/b) ' 



with F3{b) = 3sm(7r/&)(l - 9/6^)(l - 9/ Ab"^)/ sin{3TT/b). Asymptotically at large b the function 
i^3(6) ^ 1, it is really close to 1 in wide interval of argument, from 6 ~ 6 up to oo, the 
deviation does not exceed several %. However, correction due to deviation of F3 from 1 will 
be taken into account. 



Expressions (9, 10) allow to obtain the mass of multiskyrmion in simple analytical 
form as a function of parameters b and ro' 



M(B,ro,fe) 



36 sin{-K/h) 



(1 - X)iBb' + 2X)(1 - l/b') + M—^Fsib) 



(11) 



Since b ^ 2\/B, or greater, see below, we can put F3 ^ 1 and substitute at large enough B 
TT/bsin{TT/b) 1, to obtain: 



MiB,ro,b) ~ — <|5 
6n 



,AB ^ 81 



Bro + (1 - A) 



3ro 



(12) 



Now minimization over b can be done without difficulties. It provides an upper bound for 
the soliton mass, because we restricted ourselves with the profiles of the type (7), only. 
Minimizing in b is trivial, giving 



brmn = ^/ G N / G D . M {B , Tq) = — ^Ch^ 

OTT 



with 



Gn — 



Bro + {l-X) — 
3ro 



4:B 8 

Gd = ANro + (1 - A)— + AX— ^. 



(13) 



(14) 



Minimization of the mass in ro provides now r™™ and upper bound for the mass M{B). The 
extreme cases A = (SK4) and A 1 (SK6) are simple, and we present here results for both 
cases. 

The SK4 case, considered in 0. 



2 / I 
3V ^ 



1/2 



jrnin 
' ^0 



2(XM^)i/4 



M 4 /2 
"b 3^1 3 



1/2 



(2 + y/ANl/B^) 



(15) 



The SK6 case. 



81 



1/4 



limn 



M 



1/4 



(16) 
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The masses of multiskyrmions in the SK6 case at large B are smaller than those in the SK4 
case, the ratio of radii is greater: ro{SK6)/ro{SK4) ^ (6/5)1/'' ^ i.o466. 

For the SK4 variant it is possible to obtain more accurate estimates describing also 
preasymptotics in B. 

'{Bb^ + 2T}(1 - 1/62)11/2 



= 2 



and 



M{B,b)/B 



3{ANb^ + 8B) 
[{b^ + 2I/B)iANb^ + 8B)(1 - l/b^)/{3b^B)] 



(17) 



(18) 



36 sin(T:/b) 

At large enough B when it is possible to neglect the influence of slowly varying factors 
(1 - 1/62) and bsin{-K/b) we obtained P] 



6"" = 60 = 2{I/An) 



1/4 



2 / / X 



1/2 



(19) 



and, therefore, approximately 



The lower bound in (20) is taken from (4a). 

The correction to the value 60 can be found including into minimization procedure 
the factor (1 - l/fe^) and variation of bsin{n/b) ~ tt[1 - n^/{6b'^)]. It provides: 



16X3/4^]V4 

and the value b ^ bo + Sb should be inserted into (18). This improves the values of M/B 
for B = 1,2,3... but provides negligible effect for b ^ 17 and greater, since 6b ^ 1/Vb. The 
comparison of numerical calculation result and analytical upper bound (18) is presented in 
Table 1. 

For the SK6 variant from (11) we have after minimization in ro 

1/2 

min \ 2 



and 

M{B,b) 



B %sin{-K/b) 



8IF3{b)b^ ' 
5{ANb^ + 8B)_ 

{ANb'/B 



(22) 



1/4 

(23) 



which provides the upper bound for the soliton mass, for any 6, similar to (18). Combining 
with (4a), we obtain approximately 



(24) 



2IIAn\' M 8 /8AjvXF3(fe) 

3V 52 j ^ ^ 36osm(7r/6o) V 15^2 



The correction to the value of feo = ^^JB/Am can be taken into account similar to the SK4 
case, 

SbiSK,)^'-^^^^. (25) 

The results for the upper bound for the masses of the SK6 multiskyrmions calculated ac- 
cording to (23) are given in Table 2, next Section. 

A comment concerning the behaviour of the profile / at large r is necessary. It is well 
known that asymptotically at r ^ 00 / is defined by the 2~ d order term in the lagrangian 
and is proportional to / ~ l/r^ with p = 1/2 + yj2B + 1/4. So, p = 2 for B = 1, p = 3 for B = 3, 
etc. Obviously, the tail of the profile we obtained here, / l/r''/2 with b given in (15) or (16), 
is greater and falls down more slowly than the true one. This is because our purpose is to 
describe the masses of skyrmions and other global characteristics, but not the asymptotic 
behaviour of the profile. 



3 Numerical results and comparison of the SK6 and SK4 variants 
of the model 

Numerically (18) provides the upper bound for the skyrmion masses which differs from the 
masses of known RM multiskyrmions in the SK4 variant within 2%, beginning with i? = 2, 
see Table 1. The values M/BIum are calculated numerically by means of minimization of 
the functional (3). Even for B = 1, where the method evidently should not work well, we 
obtained M = 1.271 for b = 2.85 instead of precise value M = 1.232. It should be noted here 
that analytical results of paper for smallest B are slightly improved here due to better 



choice of power b{B): it is found by numerical minimization of (18), whereas in we used 
approximate formula (21) for Sb. For maximal values of B between 17 and 22 where the value 
of I is calculated, the upper bound exceeds the RM value of mass by 0.5% only. We took 
here the ratio Ri/b = in the cases where this ratio is not determined yet, the same as 
for highest B where it is known, i.e. 1.28 for SU{2) case [||], S = 32 and 64, and 1.037 for B > 6 
in ^f/(3) g. 



B 


2 


3 


4 


5 


6 


7 


13 


17 


22 


32 


64 


M/B\rm 


1.208 


1.184 


1.137 


1.147 


1.137 


1.107 


1.099 


1.091 


1.092 


1.088 


1.084 


tb 


1.45 


1.73 


1.89 


2.13 


2.30 


2.40 


3.23 


3.65 


4.15 


4.97 


6.98 


h{B) 


3.73 


4.38 


4.77 


5.35 


5.76 


6.00 


7.98 


9.01 


10.23 


12.24 


17.16 


M/B\appr 


1.227 


1.198 


1.150 


1.158 


1.147 


1.117 


1.106 


1.0976 


1.098 


1.094 


1.089 


M/B\num 


1.1791 


1.1462 


1.1201 


1.1172 


1.1079 


1.0947 


1.0834 


1.0774 


1.0766 








1.222 


1.215 


1.184 


1.164 


1.145 


1.138 


1.120 


1.115 


1.111 


1.1064 


1.101 


^B 


1.28 


1.54 


1.72 


1.88 


2.02 


2.16 


2.87 


3.25 


3.68 


4.40 


6.18 


biB)'"^'^ 


3.32 


3.96 


4.38 


4.76 


5.08 


5.43 


7.12 


8.05 


9.08 


10.86 


15.19 




1.247 


1.231 


1.198 


1.176 


1.156 


1.149 


1.127 


1.121 


1.116 


1.111 


1.106 



Table 1. The skyrmion mass per unit i3-number in universal units iir'^ Ft^ / e for the RM configurations, the 
SK4 variant, approximate and precise solutions, tb = •y/< > - mean square radius of baryon number dis- 
tribution, or the isoscalar radius, in units 2/ F-j^e. The approximate values (upper bounds) are calculated using 
formula (18) with the power b minimizing it. The muncrical values for the SU(2) model are fi'om the papers 
[D and earlier papers. The last 4 lines show the result for the SU{3) projector ansatz 1^ and approximation 
to this case, Ajy = 4/3. Calculations of Af/i?|f^'^ are made also with the help of (18) with the power b{B)^^^ 
which minimizes it. 



It is of interest to compare the same quantities (&o, fo, M/B) at large baryon numbers 
for the SK4 variant (original Skyrme model) and the SK6 variant. 

bl{SK6)/bliSK4) = [BV(XAjv)]^/2 < 1, (26) 

rUSK6)/rliSK4) = (6/5)'/' > 1. (27) 
Asymptotically at large B the ratio of upper and lower bounds 

M 4 /2\ 

Mrnin V 3 / 

i.e. the gap between upper and lower bounds is less than 4%, independently on S, the 
particular value of I and the number of flavours N. For the SK6 variant we obtain 

Rmacc/rmniSm - ^^P^\sK6 = - ( ^ ^ 1-0881, (29) 

here the gap is less than 9%, also independently on B, I and N. 
Note that the ratio 

R.ma./n.rn{SK6)/R„,a./,mn{SKA) = ro{S K6) / ro{S K A) ~ (6/5)1/^ ~ 1.0466, (30) 

so, the number (6/5)^/^^ plays a special role in comparison of both variants of the model. 

The SK6 variant reveals even weaker dependence on the quantity I than the SK4 
variant, the power 60 does not depend on I at all, cf. (15) and (16) above. 



With decreasing X the upper bounds decrease proportionally to the lower bounds. It 
should be stressed that our upper bounds for masses of the SK4 and SK6 multiskyrmions 
are obtained on a class of profile functions (7), and, probably, can be improved. 



B 


2 


3 


4 


5 


7 


13 


17 


22 


32 


64 


128 




0.9181 


0.8843 


0.8289 


0.837 


0.793 


0.781 


0.772 


0.773 


0.768 


0.763 


0.760 


m 


4.08 


4.66 


5.14 


5.55 


6.27 


8.00 


8.95 


10.01 


11.85 


16.39 


22.9 




0.951 


0.915 


0.857 


0.865 


0.817 


0.801 


0.791 


0.790 


0.784 


0.778 


0.774 


' B,RM 


1.636 


1.910 


2.047 


2.293 


2.540 


3.351 


3.767 


4.269 


5.093 


7.125 


10.037 


^^^RM I^^^R.M 


0.760 


0.747 


0.729 


0.730 


0.716 


0.711 


0.708 


0.7075 


0.706 


0.7044 


0.704 


„SK6 /„SK4 

rB 


1.128 


1.107 


1.082 


1.076 


1.056 


1.039 


1.033 


1.029 


1.025 


1.022 


1.019 



Table 2. The skyrmion mass per unit _B-nunibcr for the SK6 variant of the model (A = 1) with SU{2) flavour 
symmetry, in units iir'^Fjr/e, and comparison with the SK4 variant. The approximate upper bound M/B\appr 
is calculated using formulas (23) with the power b{B) given in third line. The isoscalar radius rB,RM in units 
2/F^e. The ratios of masses and radii tb for the SK6 and SK4 variants are presented in the last two lines. 



The masses of multiskyrmions for the SK6 variant presented in Table 2 are consid- 
erably lower than for the SK4 variant (they are between 0.76 and 0.7 of them, roughly) for 
accepted choice of parameters. They agree well with numerical results of recent paper 0. 
For the SK6 variant the B = \ configuration can be described by profile of the type (7) with 
accuracy about 3%, as for the SK4 variant: M{l)appr — 0.969 for h = 3.4 in comparison with 
numerical value Mil) = 0.940. At higher B-numbers the upper bound for masses of the SK6 
multiskyrmions is not so close to numerical values as for the SK4 variant. The difference 
is not smaller than 2%, but from practical point of view it is quite good agreement. In 
view of good quantitative agreement of analytical and numerical results the studies of basic 
properties of bubbles of matter made in and in present paper are quite reliable. 

The width (or thickness) W of the bubble shell can be estimated easily. We can define 
the half- width as a distance between points where = ±1/2, then: 

W = 'C^lni. (31) 

Looking at (6) we see that maximal value of 0' is close to 6o/(2ro), and this provides imme- 
diately W ~ 4ro/6o, in agreement with (31). Here there is some difference between the SK4 
and SK6 variants of the model, because the ratio ro/6o is different. For the original SK4 
variant at large B the thickness W ~ 2^273/^3, i.e. it is universal characteristic of all bary- 
onic bags, and does not depend on the number of fiavours N as well. For the SK6 variant 
W ~ 2{?,Ajsi1/lbB'^Y/^, i.e. it can depend slightly on B if the ratio l/B"^ has such dependence, 
and increases also with increasing number of flavours. The radius of the bubble grows with 
increasing B like ^/AnY^'^ for both SK4 and SK6 variants, see (15), (16). 

4 Toy model: the inclined step approximation 

A natural question is to what extent the "bubble" structure is a necessary property of 
multiskyrmions, and what could be instead of this. What we have is a boundary condition 
on proflle function, /(O) = tt and /(oo) = 0, and requirement to get the minimal value of the 
mass (3). In principle, the profile / could decrease according to some law different from 



(7), providing another mass and i?- number distribution, e.g. uniformly filled bag, and it 
is just the property of lagrangian (3) that the bubble structure has an advantage. A good 
illustration for this provides the toy model of "inclined step" type 0]. Let W be the width 
of the step, and ro - the radius of the skyrmion where the profile / = 7r/2. / = Tr/2- {r-ro)TT/W 
for To - W/2 < r < ro + W/2. This approximation describes the usual domain wall energy [0 
with accuracy 9.5% . 

We can write the energy in terms of W, ro, then minimize it with respect to both of 
these parameters, and find the minimal value of energy. 



on 



ttV , 2 - „ - 3X\ „y „ 3X 



This gives 

Amt^o + (1 - + 3AJ/(8r2)T 



Wmm — TT 

and the mass 



S + (l-A)3X/(8r2) 



(32) 



(33) 



Af = ^ [{Am7-1 + (1 - X)B + 3XI/{8rl)) {B + {1 - X)3I/{8rl))] (34) 

The minimization of (34) over ro should be performed now. For pure SK4 (A = O) and 
SK6 variants (A = l) it can be made trivially, and in both cases provides the same result, 
(r™")2 = ^3X/(8AAr) ~ 0.612^57^- K is close to the above result (rj,"™)2(5ii'4) ~ 0.667a/x7^v 
and (rJJ^")2(S'i^6) ~ 0.73^^/1^. In dimensional units rj,"™ = (6X/AAr)i/V(F^e). 

The thickness of the envelope W is slightly different for both models. For the SK4 
model Wrnin{SK'^ = TT [Q, I.e. it does not depend on B for any SU{N), similar to previous 
result (31) which gives W{SKA) ~ 1.8 for large B, all in units 2/(J^^e). For the SK6 variant 
Wmin{SK6) = n[3Aj^I/2B^]^/'^ > WrmniSKi) in this toy model, also similar to result obtained 
in previous Section, and increases with increasing N. Numerically, effect is not great since 
An < 2. For SU{3) group the factor is - [3An/2]^^-^ = 2^/4 ~ 1.19. 

The energy obtained in this way is 

M^.niSKA) ~ {2B + ^2,AnI/2)/Z (35) 

and 

M„MSm - f (^)''' - 0.738 5 (^)'''. (36) 

In difference from previous results, (35) and (36) do not give the upper bound for the 
skyrmion masses since some terms in expansion in W have been neglected, and for small B, 
indeed, the value of (35) is smaller than calculated masses of skyrmions. For SU{2) model 
An = I and the energy Mrmn{SK^) = {2B + ^3X/2)/3. The formula gives the numbers for 
5 = 3, 22 in agreement with calculation within RM approximation within 2 - 3% P, |^]. 

More detailed analytical calculation made in present paper confirms the results of 
such "toy model" approximation and both reproduce the picture of RM skyrmions as a two- 
phase object, a spherical bubble with profile / = tt inside and / = outside, and a thickness 
of the shell which is fixed (the SK4 model), or slightly depends on S-number and N (the 
SK6 model). 

The surface energy density can be estimated. For the SK4 model p^"''-^ ~ {2B + 
VT^)/(127rr2) ~ 3y^ANB^/I/8TT, or in ordinary units p"j^/^ - 9TTy^A^^WJlF^e/32. The average 



volume mass density in the shell is, in ordinary units, 



plf ^ ^(25 + v^)yi^F>2 ^ _.^A^ByiFy. (37) 

For SU{2) model at large B it is about (0.4 - 0.6)Gev/Fm^ depending on the value of W 
discussed above, at reasonable choice of model parameters = 0.186 Ge?;, e = 4.12 i.e. 
this density is several times greater than normal density of nuclei. For the SK6 model the 
density of matter in the shell is about ~ 0.7 - 0.8 of density for the SK4 model. 



5 The properties of the large B multiskyrmions 

It is possible to calculate analytically such characteristics of multiskyrmions as mean square 
radii of baryon number and mass distributions, tensors of inertia, etc. Let us estimate first 
the average dimensions of each cell in the envelope of the bubble of matter. The average 
area of the cell is 47rrg/(2i3 - 2), at large B it is for the SK4 model 



Scell ^^Try/T/A^/ (SB). (38a) 
Therefore, the average radius of the cell is 

rceii ^ [X/An]'^^/V3B, (39) 

i.e. it depends slightly on the ratio I/B which is close to 1 at large B. For the mentioned 
above choice of SU{2) model parameters rceii ^ 0.32_Fto, about twice smaller than radius of 
the B = 2 torus. 

For the SK6 variant the dimensions or each cell are slightly greater. 



Scell ^ iTTy/2/T5y/T/A^/B. (386) 

The radii of the baryon number and mass distributions can be calculated analytically 
for the profile (7) in terms of parameters ro, b. We have 

<r|>(ro,6) - fr'4f'dr = ^-hl-4>y/\dr = rl f (40) 
TT J ■' n J b'^cos{2-K I b) 

At large B, < r% >~ rQ[l + {2tt'^ - 16)/6^] . Evidently, when B ^ oo and 6 — > cx) then < r| >^ tq, 
as it should be expected. But for realistic values of B and b the value of < r| > is greater 
than rl. When 6 -> 4, {b^ - 16)/[b^cosi2'n:/b)] A/tt. 

Similarly we can calculate the radii of the mass density distribution. 



MiB,b)<rlj>^^ 



Ti{ANb^/4 + 2B) + (1 - X)T2{Bby2+I) + Xlb^Ts/A 



(41) 



Ti = / (1 - ^y^dr - ^^'"'^0 



b'^sin{3Tr/b) ' 

T2^ f{l- c^'rdr = „,/.7 ,.J 1 - l/b% 
J 6b'^ sin[-K I b) 

T3= /(I - (ly^)Vr^dr = .f"", fl - 1 Vl - ) . (42) 

^ J ^ ^ ' ' 15rob^sin{T:/b) \ b^ J \ Ab^ ' ^ ' 

Evidently, at large B < rf^ >^ r^, and < >~< r| >. 



It is possible also to calculate analytically the tensors of inertia of multiskyrmion 
configurations within this approximation. This will be done here for the SK4 variant which 
is now of greater practical importance than the SK6 variant. The explicit expressions for 
tensors of inertia of multiskyrmions in general SU{2) case are given in [|. They define the 
rotation energy of a skyrmion in the form 



0int 
ab 



^^a^b, 



(43) 



where angular velocities of skyrmions rotations in isotopical (wa) and usual {fib) spaces are 
defined in standard way in terms of corresponding collective coordinates and their time 
derivatives, see, e.g., and references therein. The isotopical tensor of inertia 



©ah i^ab - naUb) 



El 

4 



{df? 



+ -^diuadiub }d-^r 



(44) 



The expression for the orbital tensor of inertia is much more complicated and we shall 
not give it here, see 0] again. However, the traces of both tensors of inertia are simpler, 
especially for the RM ansatz, and depend on the quantities N = B and I H] : 



47r 



i?2 

2 



B- 



r dr. 



:|4 + ^(Br + ifi)}.^.. 



Evidently, the inequality takes place 



eL-Bel^^{j-B') J 4dr>o 



(45) 



(46) 



since I > B"^. The interference tensor of inertia is much smaller for all cases except the cases 
of spherical and axial symmetry, as for B = 2, and will not be considered here. 

The point is that at large B when multiskyrmion is close to spherical bubble the 
diagonal components of tensors of inertia can be calculated as 1/3 of corresponding traces 
(45), and off-diagonal tensors of inertia being close to zero. The accuracy of these statements 
increases with increasing baryon number. 

Now we can calculate (45) in our approximation for the profile (7). It gives: 



T2{4B + b^ 



J 

aa 



An 



[4BTi +T2(4X + B62)] 



(47) 



with Ti, T2 given above. At large baryon numbers Ti ~ 4r^/6, T2 ~ 8ro/36, and in natural units 
for tensors of inertia, 12tt'^ /FT,e^ 



0i/3 : 



So, we obtain 



(48) 



(49) 



as it should be for empty spherical bubble with its mass concentrated in its shell. The 
equality (49) does not hold, however, for the contributions of second order and Skyrme 



terms in the lagrangian separately. For the isotopical tensor of inertia we have inequahty 
6^a/3 < 2MBrl/3B. The isoscalar magnetic moment of the baryonic system is defined by the 
orbital inertia fi° ~ JB < > /3Qj, therefore we obtain from (49) 



EL 
J 



B 



(50) 



almost constant for large i?-numbers. The consideration of other electromagnetic and weak 
interaction properties of multiskyrmions is behind the framework of present paper, it will 
be made elsewhere. 

We finish this Section with a remark that the characteristics of multiskyrmions ob- 
tained here, < >, < r^/ >, Sceih moments of inertia as well as thickness and the mass 
density of the shell of the bubble given in (31) and (37) provide a complete picture of large 
B multiskyrmions as quasiclassical objects formed by the chiral fields. 

6 The role of the mass term 



Consider also the influence of the chiral symmetry breaking mass term {M.t.) which is de- 
scribed by the lagrangian 



Lm = M.t. = 



r^(l — cos f)dr, 



cosf 



(51) 



For strangeness, charm, or bottom the masses uik, tod or niB can 



be inserted instead of to^, multiplied by corresponding flavour content of the skyrmion. 
Instead of the above expression (11) we obtain now 

Mb < M{B,r^,h)^a{B,b)r^ + {\-\)P{B,h)lro + \i{B,b)lrl+mrl (52) 

with a, (3, 7 given in (11), (12) and m = 27rTO/(6sm(37r/6)). It is possible to obtain in a simple 
form the precise minimal value of the mass for the SK4 model (A = 0) 



r\ mm 

M{B, b) = — + 12TO/3 + 2a) 



at the value of 



\B,h) 



^oP- + 12 to/? — a 



6m 



1/2 



(53) 



(54) 



Eq-ns. (52) and (53) give the upper bounds for the mass of the skyrmion, because they are 
calculated for the profile (7) different from the true profile which can be obtained by explicit 
minimization of the energy functional (3) with the mass term included. In particular, it is 
well known that the tail of the profile decreases exponentially, ~ exp{-^r)/r at r > V2B/fj. 
instead of the power law. However, for multiskyrmions at large B the main contribution to 
the mass, moments of inertia of the skyrmion, etc. is due to the shell of the bubble, and the 
relative contribution of the tail (i.e. the region outside of the bubble) decreases as 1/Vb, at 
least, being not important at large B. 

When the mass m is small enough, as for the pion, the expansion in 12mi3/a'^ can be 
made, and one obtains the following reduction of the skyrmion size ro: 



3m 

''^''"-2^[a) 



I 

An 



1/4 r 



37rm 



2{2B + ^/TA^) \ A 



3/4. 



(55) 



and increase of the soliton mass 



6M = M 



m[3 
2^ 



1 - 



9m/3 
"8^ 



Mm- 



^2(25 + v^A^). 



We used that at large B 



a : 



AnB 



8B 



9tt 



Bb 
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-I ~ B' 



(56) 



(57) 



As it was expected from general grounds, dimensions of the soliton decrease with increasing 
TO. However, even for large value of m the structure of multiskyrmion at large B remains 
the same: the chiral symmetry broken phase inside the spherical wall where the main 
contribution to the mass and topological charge is concentrated Q. The value of the mass 
density inside of the bubble is defined completely by the mass term with 1-0 = 2. The 
baryon number density distribution is quite similar, with only difference that inside the 
bag it equals to zero. Eq. (52), (53) provide the upper bound for the skyrmion mass, the 
lower bound (4a) takes place again for the mass without the mass term. It follows from 
these results that RM approximated multiskyrmions cannot model real nuclei at large B, 
probably B > 12- 20, and configurations like skyrmion crystals [|lll may be more valid for 
this purpose. 

One of the issues of multiskyrmions phenomenology are the multibaryons with fiavour 
different from that of u, d quarks, in particular s, c and b fiavours, see [O, M and references 



therein. The important ingredient of the calculation of spectra of such multibaryons is calcu- 
lation of the fiavour excitation energies u)b,s, ^^b,c and The behaviour of these energies 
as a function of the baryon number is important for the conclusion if the corresponding 
flavour is more bound when B increases, or less bound. The following expression for these 
energies was obtained ( |]T2|, ^ and references therein) : 



NcB , 

^B,F = — [fJ-F,B - 1) 



86 



F,B 



with 

fiF,B ^ [1 + l6mlTBQF.B/{N,Bf]^/^ 
and the a-term and "flavour" moment of inertia Qf,b given by 



(58) 



(59) 



(1 - Cf)d\, 



eF,B = 



(1-c/) 



(60) 



Nc is the number of colours of the underlying QCD, the last term in (58), proportional to 



A^c, is due to Wess-Zumino term present in the action of the model. 



terms of relative order (F^ - F^)/mj^ because, as it was shown explicitly in 



We neglected the 



they are 



important only for small values of B, and make vanishing contribution for the large B which 
we consider here. The terms ^ (l-c/) in the integrand of (60) give contribution proportional 

^, whereas terms s} or (9F)^, etc, are due to the shell of the 



to the volume of skyrmion, 
skyrmion, only, r^. mj^ = F^mj^/F^ 
meson). When fho is sufliciently large - practically it is always fulflUed - then 



m_D is the mass of flavoured meson (k, D or B- 



For small B the energies decrease somewhat with increasing B, and this leads to the increase 
of binding energies of flavours. From the picture of large B multiskyrmions clarified here 
we can see that the quantity has the part proportional to the volume occupied by 
multiskyrmion, i.e. ^ i?^/^, whereas for Qf.b this part is considerably smaller. As a result, 
at large B the ratio 

4^ - :x (62) 



For the difference ep.B = mo - ojf.b which is important contribution into the binding energy 
of flavoured meson by SU{2) skyrmion we obtain 

where the surface, or shell, contribution to flavour inertia Q^p^^ is proportional to 1/e^ in (60). 
Only the first term in (63) comes from the volume of the multiskyrmion, but numerically it 
is negligible, except the case of the strangeness. Remaining terms are both of the surface 
(or shell) origin, in other words, the binding of heavy flavour takes place on the shell of 
the multiskyrmion. The second term, ^ Q^p^J /Tb, which dominates in magnitude, decreases 



with increasing B as and this explains the results of calculations ft O which have 



shown that the binding of flavour becomes weaker with increasing B. However, this property 
may be intrinsic for RM multiskyrmions, and can be absent for the skyrmion crystals, for 
example. 



7 Discussion and conclusions 

In and in the present paper we established the link between the topological soliton models 
in rational maps approximation, for two different modifications of the Skyrme model, the 
SK4 and SK6, and the soliton models of "domain wall" type. It was possible by means 
of consideration of a class of simple functions (7) approximating both the skyrmion profiles 
and the spherical domain walls. The ansatz (7) does not describe correctly the asymptotical 
behaviour of the profile at large distances, although provides quite accurate description 
of many characteristics of multiskyrmions. This simple picture of the RM multiskyrmions 
allows to understand some peculiarities of multiskyrmions phenomenology which appeared 
as a result of the calculations, or as computer evidence. 

The upper bound for the energy of multiskyrmions is obtained which is very close to 
the known energies of the RM multiskyrmions, especially at largest J5, and is higher than 
the known lower bound by ~ 4% only for the SK4 (pure Skyrme) variant of the model, and 
by ~ 9% for the SK6 variant. For the SK6 variant the upper bound obtained is considerably 
lower than for the SK4 variant, dimensions of solitons are larger by few %, at asymptotically 
large baryon numbers for accepted choice of model parameters. 

The following properties of bubbles of matter from RM multiskyrmions are established 
analytically, mostly independent on particular values of the quantity I: 

The dimensions of the bubble grow with B as \fB, or as X^/'', whereas the mass is 
proportional to ~ B at large B. Dimensions of the bubble decrease slightly with increasing 
N - the number of fiavours, tq - [N/{2{N -l))]'^/'^, see (15), (16). 

For the SK4 variant the thickness of the bubbles envelope (31) is constant at large B 
and does not depend on the number of fiavours, therefore, the average surface mass density 



is constant at large B, as well as average volume density of the shell (37). Both densities 
increase slightly with increasing TV. For the SK6 variant the thickness of the shell slightly 
depends on the ratio l/B"^. At the same time the mass and B-number densities of the whole 
bubble when i? —> oo, and this is in contradiction with nuclear physics data confirming 
the constant density of nuclear matter. The bubble structure of multiskyrmions is more 
pronounced for the SK4 variant of the model. The material which the shell of the bag is 
made of looks like honeycomb, or web with constant average area of each cell. 

The treatment performed in the present paper could be generalized in several direc- 
tions. More general effective lagrangians can be considered, including higher terms in chiral 
derivatives, and the bundle of profiles can be introduced instead of one, see (7): 

Pi^)^M^d, (64) 



(r/ro)^ + 1 

with the evident normalization condition / p{b)db = 1. The multiskyrmion mass can be 
written then in the form: 

M{B) = M2(S)ro + M4(S)/ro + ^^e{B)/rl + MB)/rl + ... (65) 

with fi2iB) = / p{b)a{B, b)db, etc. It would be of interest to study the advantages and prospects 
of such more general approach. Some models for higher order terms in the lagrangian have 
been considered in []T3[ for the SU{2) skyrmions. 

It follows from the above consideration that the spherical bubble or bag configuration 
can be obtained from the lagrangian written in terms of chiral degrees of freedom only, i.e. 
the Skyrme model lagrangian leads at large baryon numbers to the formation of spherical 
bubbles of matter and thus provides a field-theoretical realization of the bag-type model. In 
such models, which have been popular a time ago, the properties of the building material of 



bags have been postulated . It was the technical obstacle in formulating the bag model as 



a consistent quantum theory: the bag boundary R{t) was prescribed externally, and special 



efforts have been made to overcome this [|T5[, not successful completely. The chiral soliton 
models leading to existence of multiskyrmions are free of this drawback. It should be noted 
that there is difference of principle between hadronic bags of the MIT type and the bubbles 
of matter which appear in effective field theories as RM multiskyrmions. In the first case 
the bag is the region where the quark-gluon phase of matter is restricted, in the latter case 
it is the region with skyrmion profile function / = tt, different from its vacuum value / = 0. 
The quark-gluon phase does not enter the consideration at all, although it can exist inside 
the skyrmionic bubble, due to the known Cheshire Cat principle [1T6| , |T^ . 

This picture of the mass and i3-number distribution in the RM multiskyrmions con- 
dradicts to what is known about nuclei, however, it emphazises the role of periphery of the 
nucleus and could be an argument in favour of shell-type models of nuclei. The skyrmion 



crystals are believed to be more adequate for modelling nuclear matter. 

It would be of interest to perform the investigation of the dynamics of bubbles in the 
chiral soliton models similar to that performed recently for the simplified two-component 
sigma model, or the sine-Gordon model in (3 + 1) dimensions [r^|. Observations concerning 



the structure of large B multiskyrmions made here can be useful in view of possible cosmo- 
logical applications of Skyrme-type models, see e.g. fl^. The large scale structure of the 



mass distribution in the Universe is similar to that in topological soliton models, and 



it can be the consequence of the similarity of the laws in micro- and macroworld. We con- 
clude with a remark that analytical methods which are not typical for studies of skyrmion 
properties (see also ) allow to obtain very simple and transparent results which accuracy 
increases with increasing baryon number. 

As it was noted by referee, the parametrization of the profile function similar to (7) 
was proposed long ago in for i? = 1 hedgehog. It was f{r,ro) = 2atan{ro/r)^, or 

cosf - 



which has correct asymptotics at large r with the appropriate choice of ro, but does not 
minimize the mass of the skyrmion. Parametrization (7) with the power 6 as a parameter to 
be fitted provides description of masses especially good for the large baryon numbers. 

I would like to thank W.J.Zakrzewski, T.Ioannidou for interest in the questions dis- 
cussed in present paper and useful comments. I'm indebted to P.Sutcliffe and B.Piette for 
sending me the results of their investigations. This work is supported by the Russian Foun- 
dation for Basic Researches, grant RFBR 01-02-16615. 
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